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The electroweak 0(a) radiative corrections to the decay of the lightest MSSM Higgs boson 
to four leptons are presented, improved by the two-loop corrections provided by the program 
package FeynHiggs. We also analyze the results in the decoupling limit and investigate the 
numeric impact of contributions from the genuine supersymmetric particle spectrum. 
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1 Introduction 

Deciphering the mechanism that breaks the electroweak symmetry and generates the masses 
of fundamental particles is one of the central tasks of Tevatron and LHC. In the standard model 
(SM), the electroweak symmetry breaking is realized through the Higgs mechanism where the 
neutral component of an SU(2) complex scalar doublet acquires a non-zero vacuum expectation 
value. While in the minimal supersymmetric extension of the standard model (MSSM), two Higgs 
doublets are required, resulting in five physical Higgs bosons. Two of them, h° and H°, are CP- 
even, one is CP-odd ^4°, and the other two, H^, are charged. Among them the lightest CP-even 
Higgs boson mass is bounded from above by M h o < 135 GeV, including radiative corrections up 
to two-loop order [UIUIS]. In this mass range, the Higgs boson decays dominantly to bb pair. 
However, this is not a promising channel for the discovery of the Higgs boson at hadron colliders 
due to the large QCD background [6J. Detailed investigations of other decay modes of the Higgs 
boson are thus necessary, and such investigations also have further implications. At the LHC at 
least one MSSM Higgs boson can be discovered over all of the MSSM parameter space. In the 
region where tan (3 and M^o take on moderate values and the region with large M^o values, only 
the lightest Higgs boson would be observable [7]. In this case, precision measurements of the 
Higgs decay properties would indicate if the Higgs boson originates from the SM or the MSSM 
and for the latter case allow to derive indirect bounds on other MSSM parameters, e.g. on the 
mass of the CP-odd Higgs boson [8j. Over a large fraction of the parameter space more than 
one Higgs boson would be accessible. Then precision measurements of the lightest Higgs boson 
properties at the linear collider can help to distinguish between different soft SUSY-breaking 
scenarios El. 
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The decay properties of the SM Higgs boson have been studied intensively in the literature. 
In refs. [HlOU] the decay of the SM Higgs boson to four fermions via a gauge boson pair was 
investigated, where the complete 0(a) electroweak corrections were presented for leptonic final 
states, and for semi-leptonic and hadronic final states also complete 0(a s ) QCD corrections 
were available. The results were further improved by corrections beyond 0(a) originating from 
heavy-Higgs effects and final state radiation. 

In the present paper we consider similar processes in the CP-conserving MSSM with real param- 
eters and compute the 0(a) electroweak corrections to the decay h° — > WW* / Z Z* — > 41eptons. 
In contrast to the SM, in the real MSSM the two CP-even Higgs bosons, h° and H°, can mix 
beyond the lowest order. The resulting Higgs propagator corrections are numerically important, 
we use an effective decay amplitude to account for such corrections. In the presence of the mixing 
between Higgs bosons, the radiative corrections to the coupling of the heavy CP-even Higgs bo- 
son to gauge bosons can be numerically relevant due to heavy fermions and sfermions in the loop 
and thus should be taken into account. The numerical results are analyzed in the benchmark 
scenarios suggested in [TTJ. We also investigate the results in the decoupling limit [15] . 

This paper is organized as follows. In section [2] we describe the renormalization of the MSSM 
Higgs sector. The strategy of our computation is outlined in section [5J to section [6] and the 
numerical results are discussed in section [7] We draw our conclusions in section [S] In Appendix 1X1 
we list the analytical results for the scalar integrals that are relevant for our computation. The 
methods used to deal with the soft and collinear photon emission are briefly summarized in 
Appendix iBl 



2 Renormalization of the MSSM Higgs sector 

The Higgs sector of the MSSM consists of two complex scalar doublets with opposite hyper- 
charges, which give masses to up- and down-type fermions, respectively. The two doublets can 
be decomposed as 

H- l -[ 171 1 V2^ 1 1 " Al > \ - ^ vcos P + + i Xi) ^ 
H " = ' 772 + 7# 2 + to) J = ( vs ' m P + 7# 2 + ^ ) ' (1) 




where v = -\/vf + vj = 174 GeV and tan/3 = v 2 /vi with < j3 < ir/2. 

The tree-level masses of Higgs bosons are determined by the bilinear terms of Higgs fields in 
the tree-level scalar potential. In the <p\ and (j)2 basis, such bilinear terms give rise to a non- 
diagonal tree-level mass matrix of the neutral CP-even Higgs bosons. It can be diagonalized by 
transforming into the mass eigenstate basis h° and H°, leading to 

<°,»° - ( 1° 4 ) ■ (2) 

where the mass eigenstates h and H° are given by the following rotation 

h° 



sin a cos a \ / <pi 



H° ) \ cos a sin a 

The mixing angle a satisfies 
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with Maq the mass of the neutral CP-odd Higgs boson. This relation follows from the requirement 
that the lowest order tadpoles and the non-diagonal entries of the CP-even Higgs mass matrix 
in the h°,H° basis vanish. The mixing angles for the CP-odd and charged Higgs bosons can be 
determined analogously. The tree-level mass of the light CP-even Higgs boson, m h o , has an upper 
bound of Mz, as a consequence of the fact that Higgs self couplings are determined by gauge 
couplings in the MSSM. 

The evaluation of higher order corrections requires the renormalization of the Higgs sector. In 
the following we concentrate on the renormalization that is needed for the present work. The 
Higgs tadpole counter terms are introduced via 

T h° T h° + 5T h o , T H o -> T H o + 5T H o . (5) 

They are fixed by requiring that the renormalized tadpoles vanish, which leads to 

5T h o = -T h o , 8T H o = -T H o . (6) 

The mass counter terms <5M| and 8M^ follow from the on-shell conditions 

ReE|(Mf) = , Re£ A o(M 2 ) = , (7) 

where the superscript T denotes the transverse part of the gauge boson self energy. These 
conditions yield 

8Ml = Re S|(M|) , 8M% = Re £ A o (M%> ) . (8) 

In the mass eigenstate basis, the mass counter terms for the CP-even Higgs bosons can be written 
as 

■ ( ) ■ (9) 

Note that the mixing angles are not renormalized, these mass counter terms follow from their 
tree-level expressions as [T5] 

Sm h°H° = -2 SM z s ™2(a + P) + -6Ml sin2(a-/3) 

- [5T H o sin 3 (a - 0) - ST h o cos 3 {a-0)\ 



2Mzswcw 

Q\ 7l/f2 „„„ r,( „, i a\ i A., 

'AO 

5m 2 h0 = 5M Z sin 2 (a + 0) + 5M% cos 2 (a - 0) 



8 tan (3 cos 2 /3[M| cos 2 (a + 0) + M\ cos 2(a - 0)] , 



+ — [8T H o cos(a - 0) sin 2 (a - 0) + 8T h o sin(a - 0)(1 + cos 2 (a - 0))] 

2M z s w cw 

+ 8 tan (3 cos 2 0\M\ sin 2(a + 0) + M 2 A0 sin 2(a - 0)] , 
Sm 2 H o = 8M Z cos 2 (a + 0) + 8M\ sin 2 (a - 0) 

— —— [8T H o cos(a — 0)(1 + sin 2 (a — 0)) + 8T h o sin(a — 0) cos 2 (a — 0)] 

zMzswcw 

- 8 tan (3 cos 2 (3[M Z sin 2(q + /3) + M^ sin 2(a - /3)] , (10) 

where the counter term <5tan/3 is introduced via <5tan/3 — > tan/3 + 8t&n0, its renormalization 
condition will be given below. 

In order to have finite Green functions, the Higgs fields have to be renormalized as well. For 
the renormalization of the neutral CP-even Higgs fields, we can choose to renormalize either the 



fields h° and H or <pi and (f) 2 ■ This is in analogy to the renormalization of gauge boson fields in 
the SM, where one can renormalize either and B bosons or alternatively their mixtures, the 
7 and Z bosons. In the SM the weak mixing angle is defined by sin 2 9\y = 1 — ^w/^z m ^ ne 
on-shell scheme [20]. This defining relation is valid to all orders in perturbation theory. The weak 
mixing angle thus receives renormalization due to the renormalization of M\y and Mz- In the 
Higgs sector of the MSSM, the relations between the mixing angles and input parameters hold 
only at tree-level, the mixing angles can be kept unrenormalized. In this paper we renormalize 
the Higgs fields following [X5|. ll6j . i.e. we introduce a renormalization constant for each Higgs 
doublet, 

H x ^{\ + ^SZ Hl )H u H 2 -> (1 + ^8Z H2 )H 2 . (11) 
Their vacuum expectation values then renormalize as follows 

vx ->■ (1 + \sZ Hl ){v\ - Svt) = vi(l + \bZ Hx - — ) , 
2 2 V\ 

v 2 -> (1 + \SZ H2 ){V2 - 6v 2 ) = v 2 (l + \5Z H2 - -^) . (12) 
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The freedom of field renormalization allows us to impose the condition = The renormal 



ization of tan then follows from the renormalization of the vacuum expectation values of the 
two Higgs doublets as 



<!3) 

We can write the renormalized self energies of the CP-even Higgs bosons in terms of the unrenor- 
malized ones, the field renormalization constants and the mass counter terms 

S/,p(fc 2 ) = T, h o(k 2 ) + (sin 2 a 5Zh x + cos 2 a 8ZH 2 )(k 2 — rn^o) — Sm 2 ^ , 
S^o(fc 2 ) = T, H o(k 2 ) + (cos 2 a 5Zh x + sin 2 a 5ZH 2 )(k 2 — m 2 H o) — 5m 2 H o , 

E h o H o(k 2 ) = E h o H o(k 2 ) + sinacosa(5Zf/ 2 - 5Z Hl )(k 2 - ^(m 2 + m 2 H o)) - Sm 2 hOH0 . (14) 



The field renormalization constants 5Zh 1 and 5Zjj 2 are determined in the DR scheme. From 
Eq. (fT4^1 one finds 

*^=-[ReE / fl0 (m| n ,)| a= o] div , 

<5Z^=-[ReS; o (m 2 o)| a= o] div . (15) 



In the DR scheme 5 tan (3 is fixed by Eq. (|13p and (|15p . At one- loop level, this DR renormalization 
of tan [3 yields gauge independent results within the class of R^ gauges (the gauge dependence 
arises at two- loop level even within R^ gauges) |17| . Hence the DR scheme is a convenient choice 
for the evaluation of one-loop corrections. Moreover, this scheme has stable numerical behav- 



ior [17H19] . In this work we use the DR renormalization of tan (3 and choose the renormalization 
scale as /j, dr = 1.5M\y, which is the scale of the physical mass of the lightest CP-even Higgs 
boson for moderate tan f3 and M^o values. Other SM parameters are renormalized as in the 
on-shell scheme in |20|. 



Beyond the lowest order, the mass of the light CP-even Higgs boson receives large radiative 
corrections, predominantly from the top/stop loops, and also from bottom/sbottom loops for 
large tan f3. In addition, the two CP-even Higgs bosons can mix with each other. These lead 



4 



to finite wave function normalization factors for the external Higgs boson in our process, which 
have to be taken into account so that a correct normalization of the S-matrix is ensured. As a 
consequence of these effects, the decay amplitude of h° — > 41eptons can be written as 

M(h° -> 4Z) = ^/Z^(M h o + Z h o H oM H o) , (16) 

where the wave function normalization factors Z h o and Z h ojjo are given by [21] 



z h° 



S faO g o( fc2 ) 



l + ReS^,(fc2)-Re ^ w ^ + ^ o(fc2) , ■ -..„„ 



k 2 =M 2 n 



Z h o H o = " - v . (17) 

Myfi denotes the physical mass of h . These finite wave function normalization factors, as well as 
the physical masses of the Higgs bosons can be computed by the program package FeynHiggs [22], 
in which also the dominant two-loop corrections to Higgs boson self energies are taken into 
account. 



3 Lowest order results 

We consider the following leptonic decay processes 

h°(h) -> WW* -> e-(k 2 ) + u e {k 3 ) + fi + (h) + v^h) (18) 

and 

h ^)^ ZZ* ^e-(k 2 ) + e + (k 3 ) + Li + (h)+fi-(h) , (19) 

where the particle momenta are given in the parentheses, the helicity indices are suppressed. 
Throughout this work the masses of the final state leptons are neglected whenever possible, i.e. 
we keep them only as regulators for the collinear singularities. In these decay processes, one of 
the intermediate gauge bosons can become resonant due to the upper bound of the h° mass. The 
finite width has to be incorporated for the resonant gauge boson in order to avoid the occurrence 
of singularities. 

In this work the Feynman diagrams are generated by FeynArts [25]. FormCalc [25] and Loop- 
Tools [28j[31] are then used to algebraically simplify the amplitudes and evaluate the one-loop 
scalar and tensor integrals that do not require the finite gauge boson width as regulators. The 
computation is carried out in the 't Hooft-Feynman gauge. 



3.1 Implementation of gauge boson width 

The finite width has to be included for the intermediate gauge boson when it becomes resonant. 
However, in perturbation theory the description of resonances requires a Dyson resummation 
of self energy insertions. This mixes different perturbative orders as only partial higher order 
corrections are taken into account. As a consequence, gauge invariance might be spoiled, since 
it is preserved order by order in perturbation theory. There have been proposals for a consistent 
implementation of the width of gauge boson. One of them is the complex mass scheme |32j . in 
which the gauge boson masses are taken as complex quantities that are defined by the poles of the 
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propagators. This scheme has been successfully applied to tree-level computations as well as the 
evaluation of radiative corrections in the SM. The implementation of the complex mass scheme 
would require independent mass parameters, which is clearly not the case in the MSSM, since, for 
example, the CP-even Higgs boson masses depend on the gauge boson masses. Another proposal 
is the pole scheme [34-36J, in which one performs an expansion of the amplitude according to its 
pole structure and includes the finite width only in the pole term. The drawback of this scheme 
is that it is not applicable near and below the threshold region. In this work the width of the 
resonant gauge boson is incorporated according to the factorization scheme [37-39J, which yields 
a simple rescaling at tree-level (for convenience, we use a subscript V to denote the vector gauge 
boson, it should be replaced by W or Z accordingly for the above processes) 

k 2 - M 2 

Mb <™ = T2 M2 T V M V M bom{T V = 0) , (20) 

ky - My + iMyYy 

where ky is the four- momentum of the resonant gauge boson, M.\,orn^v = 0) denotes the decay 
amplitude before Dyson resummation. At one-loop level, the width of resonant gauge boson can 
be incorporated as follows 

ky-My ilmS^(My) 

Ml ° OP - k 2 y-M 2 v + l MyTy Ml °° p{TV ~ 0) + k 2 y-M 2 Mb ° rn 
k 2 - M 2 

— v v KA (v — m 

" k 2 _ M 2 + iM v T v J loop ' no rsX V ~ } 

^(k 2 y)-^(My) \ 

— v - ^ 2 — ^ — + OZy j M born , (21) 

where the M.i 00 p,no r.s. term in the second row denotes the one-loop corrections excluding the 
self energy corrections to the resonant gauge boson. The last term in the first row is required to 
avoid double-counting from the inclusion of finite width in the lowest order amplitude. This term 
is absorbed into the self energy corrections, yielding the last term in the above equation. In the 
factorization scheme, the non-resonant terms are treated incorrectly, since they are simply put to 
zero on the resonance. The resulting error is, however, of higher order. 



3.2 The Effective Born amplitude 

The finite wave function normalization factors induce numerically important corrections. To 
account for these corrections, we use an effective Born amplitude 

M born {Ty = 0) = (M h o(Tv = 0) + Z hOH oM° HO (Ty = 0)) 

= ^Z^M h0 (T v = 0) (1 + cot (0 - a)Z h o H o) , (22) 

where A4° (rV = 0) and A4° H0 (Fy = 0) denote the respective tree-level decay amplitude of h 
and H° before Dyson resummation. In the second row we have written the amplitude -M° H o in 
terms of -M° h o- The lowest order Feynman diagram including the Higgs propagator corrections is 
shown in Fig. [TJ 

For the process in Eq. (fT8|) . the tree- level amplitude M° h0 (Tw = 0) is given by 

2na eM\y sin(/3 — a) 1 



M u hO (T w = 0) 



s 3 k 2 - M 2 



kl-M 2 , 



[u e -(k 2 )'JpU}-v Pe (k 3 )][u K (k 5 )-/ p uj-v a+(h)} , (23) 
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Figure 1: Lowest order diagram including Higgs propagator corrections. 



while for the process in Eq. (fTUj) the tree-level amplitude M^ (Tz = 0) reads 




x [(-1 + 2sw)u e - (k 2 )^ p uj-v e + (fc 3 ) + 2s^u e - (k 2 h P uj + v e + (k 3 )] 
x [(-1 + 2sw)u „- (A; 5 )7 P u;_u M+ (fc 4 ) + 2s^u M - (k b )Yu) + v^+ (fc 4 )] 



(24) 



where we have introduced the variables fe± with fc + = &2+^3 and fe_ = A; 4 +A;5, and w-t- = ^(ldrTs). 
These decay amplitudes differ from their SM counterparts only by a factor of sin(/3 — a). 

The effective Born partial decay width is then given by 



with the squared matrix element summed over the final state polarizations and the phase space 
factor 



4 Virtual corrections 

The evaluation of virtual corrections involves several additional issues. As in the lowest order 
amplitude, the wave function normalization factors for the external Higgs boson are taken into 
account by using an effective amplitude, in which we also include the one-loop corrections to the 
coupling of H° to gauge bosons from the fermionic and sfermionic sector, since in the presence of 
mixing between the two CP-even Higgs bosons, such corrections (especially the corrections from 
the third generation fermions and sfermions) may yield sizeable contributions as they involve 
potentially large Yukawa couplings. In addition, the photonic one-loop diagrams may involve not 
only infrared singularities, but also on-shell singularities. The on-shell singularities are closely 
related to the presence of resonant gauge boson in the loop and have to be cured by including the 
finite width of gauge boson. For this purpose, the one-loop integrals that contribute to the on-shell 
singularities are computed analytically. As before, the final state fermion masses are neglected 
whenever possible in the evaluation of these integrals. If the contribution of real photon emission 
process is taken into account, the infrared singularities will cancel out. 




(25) 




(26) 
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4.1 Virtual corrections to h° — > WW* — > 4/ 



In this subsection we describe the computation of virtual corrections to h° — > WW* — > 4Z. It is 
convenient to classify the one4oop diagrams as: photonic diagrams, SM-like diagrams and genuine 
SUSY diagrams. The photonic diagrams are the same as in the SM, in Fig. [2] we show some 
examples of them. Infrared and on-shell singularities can arise only from the photonic diagrams. 
Note that the diagram with a photon exchanged between the two intermediate W bosons does 
not contribute to on-shell singularities, since only one of these W bosons can be resonant. In 
the factorization scheme, power counting tells us that only scalar integrals resulting from the 
virtual photonic diagrams can contribute to on-shell singularities. In addition, soft singularities 
occur only in these scalar integrals as well. We will evaluate these scalar integrals analytically. 
The SM-like diagrams consist of diagrams involving the SM particles other than photon and the 
MSSM Higgs bosons in the loop. Examples of these diagrams are depicted in Fig. [3l In the 
decoupling limit M^o 3> Mz, the lightest MSSM Higgs boson behaves like a SM Higgs boson and 
all other heavy Higgs bosons decouple, the contribution of these SM-like diagrams is expected 
to approach the corresponding SM contribution in this limit. The one-loop diagrams involving 
all other SUSY particles constitute the genuine SUSY diagrams, some representative of them are 
shown in Fig. [H The counter term diagrams are depicted in Fig. 

The structure of the counter term contribution from the first four diagrams in Fig. [5] has the 
same form as in the SM, while the last diagram yields (see e.g. 



M C J = M° h0 



1 5M W 5sw cos(/3 - a) 

2 sw sin(/3 — a 



5Z e + 5Z w + ^-r^f- + 1Z ^- + :zi tz f ( cos 2 /3<5tan/3 



+ - sinacosa(£Zf/ 2 — dZnS^j + 2~( s ™ 2 a ^^Hi + cos 2 a5Zn 2 

1 5M W 5sw sin(/3 — a 



M°J = Ml 



5Z e + 6Z w + -—^ + ^- + ^f- '-[ -cos^d tan/3 

2 M^r sw cos(/3 - a) 



l H° ~ JVI H° 

1 . ,„„ > \ 1 



+ - s'ma cos a(5Zn 2 — SZh^J + 2( cos2 a ^^H 1 + sin 2 a5ZH 2 ) (27) 



with 



k 2 - M 2 

The counter terms 5Zh 1 , SZh 2 and 5tan/3 have been given in Eqs. (fT5j) and (fT3|) . the remaining 
counter terms are determined in the on-shell scheme and can be found in ref. 120 . 



In the photonic diagrams, on-shell singularities can arise if the exchanged photon becomes 
soft. In the following we describe the extraction of these singularities from the virtual photonic 
contributions. 



4.1.1 On-shell singular virtual contributions 

As argued previously, only scalar integrals resulting from the photonic diagrams can contribute 
to on-shell singularities in the factorization scheme. Here we evaluate the relevant scalar integrals 
and extract from them the on-shell singular virtual contributions. These scalar integrals are 
computed for zero gauge boson width. A finite width is inserted afterwards wherever a singularity 
arises when the resonant gauge boson goes on-shell. 

For the notation of the one- loop integrals, we follow the convention of [H]. Our evaluation of 
on-shell singular scalar integrals is based on the observation that only one of the intermediate 



8 



Figure 2: Examples of photonic one- loop diagrams for the process h° — > WW* — > 41. 



gauge bosons can be resonant due to the upper bound of the mass of h°. Therefore it is possible 
to make a decomposition of these scalar integrals by extracting from them the non-resonant gauge 
boson propagator. This decomposition leaves the on-shell singular parts of the integrals intact 
but simplifies the extraction of such singular parts from the integrals. 

We start from the on-shell singular scalar 5-point integral. Assuming the gauge boson with 
momentum is resonant, the scalar 5-point integral arising from the photonic pentagon diagram 
in Fig. [2] can be decomposed as follows 

Eo(-k 2 , -k+, k-,k 4 , A, m e , M v ,M v ,m^) 

D (-k 2 , k- , k 4 , A, m e , My ,m M ) - D (k 2 + k 4 ,-k 3 ,k 2 + k^,m e ,m^,M v ,M v ) 

- S 23 E l (-k 2 ,-k + , k-, fc 4) 0, m e , My, My, m M ) 

+ (S 24 + S u )E 2 (-k 2 , -k + ,k-,k 4 ,0,m e ,M v ,M v ,m^) 

- 2S 23 E 3 (-k 2 , -k + , k- , fct, 0, m e , M v ,M v ,m^) 

+ (S 24 + S 25 + S u + S 35 )E 4 (-k 2 ,-k + ,k„,k 4 ,0,m e ,M v ,M v ,m fl )} , (29) 

where Sij = (h + kj) 2 , and a fictitious photon mass A is introduced to regularize the soft singu- 
larity. While the second scalar 4-point integral in the curly bracket is finite, the first one contains 
both the on-shell singular logarithm ln(/c^ — My + if) and the resonant factor l/(k 2 _ — My), 
as well as the soft singularities, as one can see from its analytical expression in Appendix |Aj 
The tensor coefficients E{ that arise from the covariant decomposition of the vector 5-point in- 
tegral (see ref. |41j ) do not contain soft singularity. They may contain the singular logarithm 
ln(fe^. — M v + ie), but not the resonant factor l/(k 2 _ — My). Therefore it follows from Eq. (|2ip 
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Figure 3: Examples of one-loop diagrams for the process h° — > WW* — > 41 involving SM particles 
other than photon and the MSSM Higgs bosons. 

that only the first scalar 4-point integral gives rise to a singular logarithm ln(A^ — My + ie) when 
approaches My. This singular logarithm has to be replaced by ln(/c^ — My + iMyTy) in the 
final result, as will be done later on for other scalar integrals involving this singular logarithm. 
This replacement does not disturb gauge invariance [36|.l43j . 

The photonic box diagrams in Fig. [2] involve scalar 4-point integrals that might contribute to on- 
shell singularities. These integrals can be decomposed analogously by extracting the non-resonant 
gauge boson propagator, yielding 

D (-k 4 , k+, , 0, m M , My, My) 

C (-k 4 , -k-, 0, m M , My) - C (-k 5 , fc 4 + A;+, m M , My, My) 

+ (<S 24 + Ss^D^-h, k+, -fc_, 0, m MJ My, My) 
+ (S 2 4 + S u + S 25 + S 35 )D 2 (-k 4 , k+, -k- , 0, rrif^My, My) 

-2S 2 3D 3 (-h,k + ,-k-,0,m^Mv,M v )} , 
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Figure 4: Some representative genuine SUSY one- loop diagrams for the process h° — >• WW* — > 41. 




Figure 5: Counter term diagrams for the process h° — > WW* — > 41. 



D {-k 2 , fc_, -k+,0, m e , My, My) 

= u2 ^ M 2 { C o(-k2, fc_,0,m e ,Mv) - C (-&3, fc 2 + L,m e ,M^,My) 

- S 23j Di(-/c 2 , fc -> _A; +' °> m ^F.^) - 2S 23 D 2 (-k 2 , fe-, -A;+, 0, m e , M v , My) 

+ (S 2 4 + S M + S 25 + S 35 )D 3 (-k 2 ,k^,-k + ,0,m e ,M v ,Mv)} . (30) 



The analytical expressions for the on-shell singular scalar 3-point functions in Eq. (|30p can be 
found in Appendix [XI 

With these decompositions, the on-shell singular virtual contributions can be extracted straight- 
forwardly. As only scalar integrals contain soft singularities, the analytical results of the scalar 
integrals can also be used to extract soft singularities from the virtual contributions and allow 
an analytical check of their cancellation when combining with the real corrections. The soft and 
on-shell singular terms arising from the photonic box and pentagon diagrams can be summarized 
in a correction factor to the tree-level amplitude 



* a sing KA Q rsing . M ( a f (S 2 4 + S 34 )(k_ - My) 

M h %, p -M h0 S bp » = Mio{-—{ k 2 + -M v C {-h,-k_,0,m u ,M v ) 

+ (S 24 + S 25 )C (-k 2 , k-, 0, m e , My) 

+ S u {k 2 _ - M v )D (-k2,k-,h,\,m e ,My,m^} . (31) 

The photonic vertex and self energy diagrams also contribute to soft and on-shell singularities. 
These contributions originate from the photonic corrections to the Wff vertex, the W boson 
self energy and the field renormalization constants of the external charged fermions. The field 
renormalization constants of W boson give rise to soft singular contributions as well. However, 
they only appear in intermediate stages and cancel out in the full matrix element. Owing to the 
presence of one additional resonant propagator, the photonic W boson self energy insertion gives 
rise to a correction factor to the lowest order amplitude of the form Bo(kC., 0, My)/{k_ — My). 
After the W boson mass renormalization, the correction factor is modified to be proportional to 



B (k 2 _,0,M v )-B (M v ,0,My) _ 1 / k 



ie) , (32) 



k 2 _ -My k 2 V M y 

which is on-shell singular and has to be regularized by the width of W boson. The photonic 
vertex correction gives rise to a correction factor involving the on-shell singular scalar 3-point 
integral. Putting all these together, one finds the following correction factor that contains the 
soft and on-shell singularities from the photonic vertex and self energy diagrams 

M h o™<> =M° h0 5 s ™ 9 = M° h0 {^(klC (-k 2 ,-k + ,0,m e ,M v ) 



£<*-**/*Wi-£ 



«)} 



(5-M£/Ad)ln( 1-^73-- ie) \ > 



where the subscript IR denotes the infrared singular part of the field renormalization costants, 
whose expression can be found, e.g. in [20] . The correction factors defined in Eqs. (|3ip and 
(|33p include all the soft singularities from virtual photonic diagrams. The collinear singularities, 
however, are not fully contained in these correction factors, since the tensor integrals that are not 
accounted for in these factors, contain collinear singularities as well. 
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In the case that the other gauge boson becomes resonant, the correction factor resulting from 
the box and pentagon diagrams becomes 

+ (S24 + S 3 4)C (-h, k+, 0, rrifi, M v ) 

+ S 2i (kl - M v )D (-k 2 , -k+, h, A, m e , My , m„)) } . (34) 

The correction factor from the self energy and vertex diagrams can be obtained from Eq. (|33p by 
replacing k 2 _ with k\ in the last term. 

Now we are able to extract the on-shell singular virtual contributions from Eq. ()3ip . (|33p and 
the results of the scalar integrals in Appendix^] This gives rise to the following correction factor 

- <4i - [f (^) V § - 1 } - f - *) • <*> 

In the case that the other gauge boson becomes resonant, the correction factor describing on-shell 
singular virtual contributions can be obtained by making the following replacement in the above 
equation 

S25 o S 34 , k 2 _ k\ . (36) 



4.1.2 Soft and collinear singular virtual contributions 

For the process h° — > WW* — > Al, the photonic diagrams do not build a gauge invariant subset 
by themselves and their contributions are UV divergent. One can extract the soft and collinear 
singularities and on-shell logarithms from them and combine with the real photon bremsstrahlung 
to build the QED-like corrections, which are both IR and UV finite. The soft and collinear 
singular contributions from the virtual corrections can be extracted by making use of the well- 
known Kinoshita-Lee-Nauenberg (KLN) theorem [44], according to which the soft and collinear 
singularities are canceled out between the real and virtual corrections for sufficiently inclusive 
quantities. In the decay of h° to leptons, there is no initial state radiation of photons. The soft 
and collinear singular parts of the virtual corrections are exactly given by the singularities in the 
final state photon bremsstrahlung, but with opposite sign. The latter can be computed, e.g. with 
the dipole subtraction approach [46 , 48H5U] . Following [SD], the soft and collinear singular parts 
of the virtual contributions can be defined as 

5 5 2^2 
dTvirt,sing = dT born ^- ^ ^ Qi Qj \L(Sij,m%) + L{Sij,m]) ^- + 3) , (37) 

i=2j=i+l 

where dTi, orn is the lowest order decay width, Q, L denote the charge of final state fermions, and 
the function L(Sij,m^) is given by 

Note that when computing these contributions, we do not include the soft and collinear singular 
parts arising from the virtual photonic corrections to the mixed tree-level amplitude, i.e. to the 
second term in Eq. (|22p. since the corresponding photonic loop diagrams are not included in the 
virtual contribution. The definition of the soft and collinear singular virtual contribution is, of 
course, not unique. An alternative definition can be found, e.g. in [5T] . 
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The IR and on-shell singular virtual contributions combined with the contribution of real photon 
bremsstrahlung yield the QED-like correction. As mentioned, the soft singularities from photonic 
virtual diagrams are fully contained in the correction factors 5^ and 5^™ defined in Eqs. (|3ip 
and (|33p . Subtracting from them the virtual singular factor in Eq. (|37p (note that one has to 
take twice the real part of the correction factor in (|3ip and (|33p ). the remnant must be free 
of soft singularities. This provides an analytic check on the cancellation of soft singularities. 
The cancellation of collinear singularities is checked numerically, since the collinear singularities 
appear in scalar integrals as well as in tensor ones, the contribution of the latter is not accounted 
for in the correction factors. 



4.2 Virtual corrections to h° -> ZZ* ->■ 4/ 

The computation of virtual corrections to the decay of hP to four leptons via a Z boson pair can 
be carried out analogously. Consider the process 

h^h) -+ ZZ* -> e-(k 2 ) + e + (k 3 ) + fi + (k 4 ) + n~(k 5 ) . (39) 

We also define fe+ = k 2 + k% and = + k§ as before. The virtual one-loop diagrams can be 
classified as for the process h — > WW* — > 4/, and the h°ZZ/H°ZZ counter term contributions 

can be easily obtained from Eq. ([27]) by simple replacements: 5Z W -> 5Z Z , ^ ^(l - 2^f) , 
while the structure of the other counter terms is as in the SM. In Fig. [6]we show only the photonic 
diagrams that contribute to the on-shell singularities. 

In this decay process, the intermediate gauge bosons are neutral, thus only pentagon diagrams 
with a photon exchanged between two external charged fermions can contribute to on-shell sin- 
gularities. The scalar 5-point integrals can again be decomposed by separating the non-resonant 
gauge boson propagator. Here we give explicitly the correction factor resulting from the photonic 
pentagon diagrams that contains the soft and on-shell singularities (assuming the gauge boson 
with four-momentum k- is at resonance) 



csing a ( 

p 2 A 



Li 2 



S 25 + k 2 _-M 2 



'24 



V m e m^ J \ XMy 

- In 2 ffl - In 2 ( M v-k 2 --S 2 ,-S 2 , _ x _ ^ 
\M V J V m e M v J 3 

+ (2 v> 3, 4 o 5) - (2 o 3) - (4 <-> 5)} , (40) 

where the expressions in the parentheses are obtained from that in the squared bracket by inter- 
change of indices. In the case that the other gauge boson is resonant, the correction factor can 
be obtained from Eq. ()40p by making the following substitutions 

S 25 ^S 3 4, k 2 _^k\, rrie^m^. (41) 

The on-shell singular virtual contributions can be easily extracted from these correction factors. 
The photonic corrections to the Zff vertex and the field renormalization constants of the external 
charged fermions contribute to soft singularities, but not to on-shell singularities. There is no 
contribution to these singularities from Z boson self energy corrections. 
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Figure 6: Photonic one-loop diagrams that contribute to the on-shell singularities. 
4.3 Application to semileptonic and hadronic final states 

In the discussions above, we only consider the decay of h° to leptonic final states. The procedure 
outlined there can in principle be applied to semileptonic and hadronic final states as well. For 
semileptonic and hadronic final states, the number of photonic diagrams involving on-shell singu- 
larities would increase, but their evaluation can be carried out in complete analogy to what we did 
for leptonic final states. With appropriate substitution of momenta and masses, the analytical 
results for the scalar integrals in Appendix [A] can still be used to extract the on-shell singularities 
from the virtual contributions. The soft and collinear singular virtual contributions are again 
given by Eq. (f3T|) . The only possible exception in which the procedure used for our computation 
can not be straightforwardly applied is that the final state involves heavy down-type fermions, 
e.g. b quarks. In our computation we keep the final state fermion masses only as regulators for 
collinear singularities and neglect them elsewhere. If the final state involves b quarks, we can not 
neglect their masses, since their coupling to Higgs bosons can be enhanced at large tan (3 and give 
rise to numerically important effects. 



5 Real corrections 

In order to achieve an IR-finite physical result, the combination of virtual and real corrections is 
required, as a consequence of the fact that the experimental resolution of soft photons is limited. 
When considering the corresponding real photon emission for the processes given in Eqs. (|18[) 
and (|19p . the final state fermion masses are consistently neglected unless they have to be kept as 
regulators. Soft and collinear real photon emission give rise to singularities, which are regularized 
by the photon mass A and the fermion mass mj, respectively. These soft and collinear singularities 
are treated within two different approaches, i.e. the phase space slicing and the dipole subtraction 
approaches. They will be briefly summarized in Appendix iBl 

When evaluating the hard Bremsstrahlung contribution a constant width is introduced for each 
gauge boson propagator. For the propagator that is not resonant, the corresponding error is of 
higher order and negligible. If the collinearly emitted photon is not treated inclusively in the 
evaluation of distributions, the logarithm involving the light fermion masses would survive and 
become visible. 
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6 Higher order final state radiation 



The emission of photons collinear to final state charged fermions leads to corrections that are 
enhanced by large logarithm involving the fermion masses. If the collinear ly emitted photon is 
treated inclusively, this logarithm will cancel out as a consequence of the KLN theorem. If this is 
not the case, for instance, in the evaluation of distributions of final state muons, this logarithm will 
survive and might yield sizeable effects. Therefore one should take into account the corresponding 
higher order corrections. This can be done by the structure- function method [52, 53] based on 
the mass factorization theorem. Here we incorporate the effects of the higher order final state 
radiation following [9] [50] and choose the factorization scale as the physical mass of the light 
CP-even Higgs boson mass. 



7 Numerical results 

For the numerical evaluation, we use the following inputs for the SM parameters |58y59j 

Gf, = 1.16637 x 1(T 5 GeV~ 2 , a = 1/137.03599968 , 

M w = 80.403 GeV , T w = 2.141 GeV , 

M z = 91.1876 GeV , T z = 2.4952 GeV , 

mt = 172.7 GeV , m b = 4.2 GeV . 

The lowest order matrix element is parametrized in such a way that it absorbs the running of the 
electromagnetic coupling and the leading universal corrections to the p parameter, i.e. we use 
the effective coupling derived from the Fermi constant 

= V2G„Mfc* w 

TT 

for the Born amplitude. In the relative 0{a) corrections, we use the coupling a = a(0). To avoid 
double-counting from using a G ^, the charge renormalization constant in Eq. (|27p is modified to 

5Z e = 5Z e - ^Ar , (43) 

where Ar summarizes the radiative corrections to the muon decay. In the evaluation of distri- 
butions, a real photon closer than 5° to the emitting charged fermion or with energy less than 
1 GeV is combined with the emitting charged fermion in the inclusive treatment. 

In this analysis we investigate the results in several suggested benchmark scenarios [llj . which 
are defined so that the two parameters that govern the tree-level Higgs sector, M^o and tan/3, are 
varied while the other parameters that enter via radiative corrections are fixed. In these scenarios 
a common soft SUSY-breaking parameter M SUSY , as well as the same trilinear coupling for the 
third generation slepton and squark, is chosen for simplicity. The U(l) gaugino mass parameter 

is given by the GUT relation M\ = \-^-Mi- Throughout the parameter scan, the experimental 
mass exclusion limits from direct search of super symmetric particles and the upper bound on 
the SUSY corrections to the electroweak p parameter [58J have been taken into account. In the 
parameter region 50 GeV < M^o < 1 TeV and 1 < tan (3 < 50, the bound derived from the 
BR(i? — > X s ^y) prediction has ruled out the gluophobic scenario [60], therefore we will not discuss 
this scenario here. The investigated scenarios are listed as follows: 
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1. The m™ ax scenario 

The parameters in this scenario are given by 

M SUSY = 1 TeV , fj, = 200 GeV , M 2 = 200 GeV , 

X t = 2M SUSY , A b = A t = A T , m- g = 0.8M SUSY . (44) 

where Xt is the mixing parameter of the top squark sector and m g is the gluino mass. This 
scenario yields a maximal value of the lightest CP-even Higgs boson for given M A o and tan/3. 

2. The no-mixing scenario 

The only difference of this scenario from the m™ ax scenario is the vanishing mixing in the top 
squark sector and a higher value of M SUSY , where the latter is chosen to avoid the exclusion 
bounds from the LEP Higgs searches [61|,I62J . The parameters in this scenario read 

M SUSY = 2 TeV , n = 200 GeV , M 2 = 200 GeV , 

X t = 0, A b = A t = A T , m g = 0.8M SUSY . (45) 

3. The small-a e ff scenario 

In this scenario a suppression of the hPbb coupling can occur. The parameters are given by 

M SUSY = 800 GeV , y. = 2.5M SUSY , M 2 = 500 GeV , 

X t = -1100 GeV , A b = A t = A T , m g = 500 GeV . (46) 

In the decay processes considered in the present work, only SM particles are involved in the 
final state. It is interesting to compare the SM and the MSSM predictions for the partial decay 
widths. As discussed previously, in the limit that the mass parameter M A o is much larger than 
the electroweak scale, all the heavy Higgs bosons will decouple and the contribution to the partial 
decay width from loop diagrams involving the SM particles and Higgs bosons will approach the 
SM prediction for a Higgs boson with the same mass. In order to compare the partial decay 
width of hP to four leptons in this limiting case with the SM result [3], we also carry out the 
computation with the input parameters of ref. [9], and choose the SUSY parameters M SUSY , \x 
and M 2 to be M SUSY = \i = M 2 = M A o, so that the supersymmetric particles decouple when 
Mjifi becomes large. The remaining SUSY parameters are chosen as in the m™ ax scenario. In 
Fig. [7] we show the one-loop corrected partial decay width of h to four leptons excluding the 
contribution of genuine SUSY diagrams as a function of M A o . In the limiting case that M A o gets 
large (M40 > 1.5 TeV), we find an agreement with the SM results. 

If the generic mass scale of SUSY particles M SUSY is not much larger than the electroweak scale, 
the genuine SUSY particles will not decouple even in the limit M A o 3> Mz- To investigate the 
numeric impact of contributions from the genuine SUSY particle spectrum, we compare the one- 
loop corrected partial decay width of hP with and without the genuine SUSY loop contributions. 
In Fig. [8] the lowest order and corrected partial decay widths of hP to four leptons are depicted for 
-^susy = jU = M2 = M40, with the remaining parameters chosen as in the m™ ax scenario. The 
numerically most important one-loop corrections have been incorporated into the lowest order 
result by using the effective amplitude and the effective coupling. For the decay hP — > e~i> e y + v^ L , 
the relative size of the full MSSM loop contribution varies between —2.5% and —2% for tan /3 = 5, 
while it varies between —2% and —1.5% for tan/3 = 29. For the process hP — > e~e + /i + /i~, the 
relative corrections change from —3% to —2% and from —4% to —3% for tan /3 = 5 and 29, 
respectively. As can be seen from the plots, for both processes, at large M A o the blue curve that 
includes the genuine SUSY loop contributions and the red curve that does not almost coincide 
with each other, indicating that the effects of the genuine SUSY loop contributions are negligible 
in the decoupling limit. In Fig. Owe choose a relatively small value 300 GeV for M SUSY , so that the 



17 




500 1000 1500 2000 2500 500 1000 1500 2000 2500 

Figure 7: Corrected partial decay width of h° from diagrams excluding genuine SUSY particles 
as a function of M^o , with M SUSY = [i = M2 = M^o and tan j3 = 6. The remaining parameters 
are chosen as in the mP^ scenario. 
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Figure 8: Partial decay width of h° as a function of M^o, with M SUS y = M = -^2 = M^o and 
tan f3 = 5, 29. The remaining parameters are chosen as in the m^ iax scenario. The dashed line 
denotes the tree- level result. The blue line shows the corrected width with full MSSM corrections, 
the red line shows the corrected width excluding genuine SUSY contributions. 
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Figure 9: Partial decay width of h° as a function of M A o, with M SUSY = 300GeV and tan /3 = 5, 29. 
The remaining parameters are chosen as in the m™ ax scenario. The right plot shows the effects 
of the genuine SUSY contributions. 



genuine SUSY spectrum is not too heavy. The relative corrections vary from —2.5% to —1.5% 
and from —3% to —1% for tan/3 = 5 and 29, respectively. However, the genuine SUSY loop 
contributions again yield negligible effects in the large M A o limit (while for small value of M A o, 
their contributions can reach several percent), as one can see from the right plot of Fig. [9j This 
implies that for the processes investigated here the decoupling behavior is essentially dominated 
by the mass parameter M A o, thus it will be rather difficult, even if one-loop corrections are taken 
into account, to distinguish the lightest MSSM Higgs boson from the SM Higgs boson if M A o is 
large. 

We also perform a generic scan over the most relevant parameters of the Higgs sector, M A o and 
tan /3. Fig. [101 shows the results for the one-loop corrected partial decay width of h° — > e~D e [i + v^ 
including the full MSSM corrections in three different scenarios, where the corrections to the 
H°WW vertex from loops involving heavy fermions and sfermions are not included. For M A o > 
500 GeV, the results hardly vary with M A o and therefore are not shown there. For small M A o 
values (M A o < 140 GeV), the decay width is rather small due to a cancellation between the two 
parts of the effective Born amplitude Eq. (|22p . When M A o and tan/3 increase, the Higgs boson 
mass and thus the decay width increases rapidly in all three scenarios and reach a plateau after 
tan/3 > 15 and M A o > 220 GeV. In the small-a e g scenario there is a slight decrease with M A o 
for moderate and large tan/3 and M A o > 220 GeV. This is basically due to the slight decrease 
of the light CP-even Higgs boson mass with M A o in this region. The relative corrections are 
— 1.5% ~ 1% and —3% ~ —2% respectively in the and no-mixing scenarios; in the small- 

er scenario they do not exceed —4% unless for small M A o values (M A o < 140 GeV), where the 
cancellation between the two parts of the effective Born amplitude Eq. (|22p can yield a rather 
small lowest order result, and the size of the radiative correction is comparable to the lowest 
order result. This is not shown in Fig. [10] so that the generic size of the relative corrections can 
be clearly seen. 

Owing to the mixing between the two CP-even Higgs bosons, the coupling of H° to gauge bosons 
needs to be taken into account as well. Although at tree-level this coupling is usually suppressed, 
the one-loop contribution can be numerically relevant, since the fermionic and sfermionic loops 
involve potentially large Yukawa couplings. In Fig. [TTJwe show the correction to the partial decay 
width of h° — > e~v e ^x + Va due to the third generation fermionic and sfermionic loop contribution to 
the H°WW coupling. Such correction involves both the Yukawa couplings and the wave function 
normalization factors resulting from the mixing between the two CP-even Higgs bosons. While 
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the latter can lead to a suppression/enhancement to the correction when M^o/tan/? increase, 
the former tend to enhance the correction when tan /3 increases. The combination of these effects 
may lead to local extremum in the M A o — tan /3 plane. For our scan, a maximum appears at 
tan P = 9, M A o = 140 GeV in the mjf ax scenario and at tan /3 = 19, M A o = 120 GeV in the no- 
mixing scenario. In the small-a e g- scenario, there is a maximum at tan/3 = 9, M A o = 140 GeV. 
The relative corrections are positive and stay below 1.5% in the m™ ax and no-mixing scenarios, 
and vary from —0.5% to 1% in the small-a e fj scenario. In all three scenarios the contribution of 
these fermionic and sfermionic loops decreases very rapidly when M A o becomes large. 

In Fig. [T2] and [13] we show the same plots for h° —> e~e + fi + fi~ , which exhibit similar features 
as Fig. [TU]and[TT] As shown in Fig. \T2[ the relative corrections are negative in all three scenarios. 
In the m™ ax scenario the relative correction varies from —3% to —0.5%, while it stays between 
—5.5% and —3% in the other two scenarios. In the small-a e fj scenario where the cancellation in 
the effective Born amplitude can occur, the relative correction can reach 95%. We also find that 
there is a maximum at tan /3 = 11, M A o = 140 GeV in the m^ iax scenario. The maximum occurs 
at tan/3 = 23, M A o = 120 GeV in the no-mixing scenario and at tan/3 = 9, M A o = 140 GeV in 
the small-a e g- scenario. 

Fig. Q31 shows the invariant mass distribution of the /i + z^t pair in the decay process 
e~v e \i r v ii in the m™ ax scenario, where the parameters are chosen as tan/3 = 30, M A o = 120 GeV 
and 400 GeV. As the mass of h° stays below the production threshold of the gauge boson pair, 
only one intermediate gauge boson can be resonant. From the plots it can be seen that in addition 
to the peak around the W boson mass, there is another broad peak at a small invariant mass. 
This is the point where the other W boson gets resonant. In the right plot the broad peak is 
closer to the W resonance peak, as the Higgs boson mass is larger. In Fig. [15] we show the 
relative corrections to the invariant mass distribution of the pair in the m™ ax scenario. 

From the plots one can find an enhancement at low invariant mass due to the emission of photon 
off the final state fermion. In the case that the collinearly emitted photon is not combined with 
the emitting fermion, the logarithm involving the light fermion mass would survive and give rise 
to large corrections. This is shown by the blue curves in Fig. PT5| where the red curves show 
the results with photon combinations, i.e. the collinear photon is combined with the emitting 
fermion. If the fermion masses are consistently neglected, the invariant mass distribution of the 
e~v e pair with collinear photon combination will coincide with the red curves in the plots. In 
Fig. US we show the relative contributions due to the higher order final state radiation, and due 
to the corrections of the third generation fermions and sfermions to the H°WW coupling, where 
in the right plot the latter is not shown since it is strongly suppressed by the wave function 
normalization factors and is completely negligible for large M A o values. For both M A o values, 
the higher order final state radiation can lead to corrections less than 2%. The corrections of 
the third generation fermions and sfermions to the H°WW coupling give rise to a contribution 
less than 1% for M A o = 120 GeV. In Fig. [T7]the invariant mass distributions of the /J ,Va pair 
without photon combination in the no-mixing and small-a e g- scenarios are shown with tan /3 = 30 
and M A o = 400 GeV. The relative corrections can reach ~ 25% in both scenarios. 

Fig.[T5]to[22show the corresponding invariant mass distributions for the process h° — > e~e + n + /j,~ 
in three benchmark scenarios. In the left plot of Fig. \18\ the broad peak at low invariant mass is 
not clearly visible in the depicted region because of the larger mass of Z boson compared to the W 
boson mass. As shown in Fig. (20] for both M A o values, the higher order final state radiation gives 
rise to larger corrections than in the previous process, as the final state now involves two muons. 
While the correction to the H°ZZ coupling leads to a contribution of ~ 1% for M A o = 120 GeV 
in the m™ ax scenario, it is completely negligible for M A o = 400 GeV. In Fig. [21] the invariant 
mass distributions in the no- mixing and small-a e ff scenarios are depicted with tan/3 = 30 and 
M A o = 400 GeV (no photon combination). In both scenarios, the relative corrections vary from 



20 



30% to ~ 10%. 



8 Conclusions 

We have investigated the leptonic decay processes of the light CP-even Higgs boson in the MSSM 
via a gauge boson pair and computed the corresponding 0(a) electroweak corrections, improved 
by the two-loop corrections provided by FeynHiggs. 

The tree-level coupling of the lightest CP-even MSSM Higgs boson to the SM fermions and 
gauge bosons approaches that of the SM Higgs boson in the decoupling limit M^o 3> Mz- At 
one-loop level, the radiative corrections that arise from loops involving the SM particles and 
the MSSM Higgs bosons also tend to the SM result in the decoupling limit, since all heavy 
Higgs bosons decouple in this limit. We computed such corrections in this limit and found an 
agreement with the SM result. The decoupling of the genuine SUSY particles is governed by 
their characteristic mass scale. If this mass scale is not much larger than the electroweak scale, 
these genuine SUSY particles might yield sizeable effects. However, our results show that they 
only yield negligible effects in the limit that M^o 3> Mz, even for a relatively light genuine SUSY 
spectrum. This indicates that in our case the decoupling behavior is essentially dominated by the 
mass of the CP-odd Higgs boson mass, thus it will be rather difficult, even if one-loop corrections 
are taken into account, to distinguish the lightest MSSM Higgs boson from the SM Higgs boson 
if M A o is large. 

For the decay processes considered in this work, the relative corrections to the decay width 
turn out to be of the order of several percent in the investigated SUSY scenarios unless for small 
CP-odd Higgs boson masses in the small-a e fj scenario, where a cancellation in the effective Born 
amplitude can occur and leads to a small lowest order result and a large relative correction. The 
corrections to the distributions are significant. Owing to the mixing between Higgs bosons, the 
potentially large corrections from the third generation fermions and sfermions to the coupling of 
the heavy CP-even Higgs boson and gauge bosons have been taken into account as well. Such 
corrections are comparable to other one-loop corrections in all investigated scenarios at small CP- 
odd Higgs boson masses, while decrease quite rapidly and become completely negligible when the 
CP-odd Higgs boson mass becomes large. In the evaluation of distributions, we also discussed 
the contribution of the higher order final state radiation, the effects of which are typically of 
the order of several percent and comparable to the contribution from the heavy CP-even Higgs 
boson-gauge bosons vertex correction at small values of the CP-odd Higgs boson mass. 

The analytical expressions of the on-shell singular scalar integrals listed in Appendix [X] can be 
useful for the computation of other processes involving unstable particles in the loop. 
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A Results of integrals 

In this appendix we list the analytical expressions of the scalar integrals involving on-shell singu- 
larities discussed in section 14.1.11 The variables ki , Sij and the mass parameters are also defined 
there. These integrals are obtained with the help of of refs. |631 f64j. The expressions of inte- 
grals are given for zero gauge boson width. The real squared mass My should be replaced by 



21 
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where A = — 7_e + ln(4-7r) and /i is the reference mass scale of dimensional regularization. 



B Treatment of real soft and collinear photon emission 

In this appendix we briefly describe the two approaches used to deal with the soft and collinear 
photon emission: the phase space slicing and the dipole subtraction method [481150] . In both 
approaches we regularize the soft and collinear singularities by the photon and fermion masses, 
respectively. 
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B.l Phase space slicing 



For the real photon Bremsstrahlung the phase space integral diverges in certain regions. One 
can divide the phase space into singular and non-singular regions. In the non-singular region the 
integral is finite and can be evaluated numerically without regulators. In the singular region the 
integral has to be evaluated analytically with regulators. The singular region consists of the soft 
region, where the photon energy is smaller than a given cutoff AE; and the collinear region, in 
which the photon is emitted collinearly (but not soft) to a charged fermion, namely the angle 
between the emitted photon and the charged fermion is smaller than an angular cutoff AO. The 
real corrections can be decomposed as follows 



d T^ = dT so f t + dT coll + dY h ^ 



(48) 



In the soft and collinear regions, the squared matrix element |7\4^°^ 4 ' 7 | 2 factorizes into the lowest 
order squared matrix element and a universal soft or collinear factor. The five particle phase space 
also factorizes into four particle phase space and a photon part, so that the photon momentum can 
be integrated over analytically. In the soft region, the soft photon approximation can be applied, 
in which the photon 4-momentum is omitted everywhere except in the IR singular propagators. 
Note that the fermion masses are kept only as regulators for the collinear singularities, the soft 
photon correction factor can then be written as [50] 
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(49) 



In the collinear region, the squared matrix element and the phase space also factorize as in the 
soft region, and the collinear factor that describes the collinear final state radiation is given by 



dT co ii = dTbo 



y±4 
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3 , (AE 
2 + \~W 



V m% J 




+ 3 - — > ■ (50) 



The cutoff parameters AE and AO should be chosen sufficiently small so that the soft photon 
and leading-pole approximations apply. On the other hand, they should not be too small so that 
the instabilities of numerical integration can be avoided. Also note that this result assumes that 
a photon emitted collinearly to a charged fermion is treated inclusively, namely it is combined 
to the emitting charged fermion. As a result, all dependence on the photon and fermion masses 
will drop out in the final result. If the collinear photons are not treated inclusively, then in the 
collinear region one has [65] 
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(51) 



with Zi = k®/k® and the splitting function 



l + z 



(52) 
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Here kf and k® denote the energy of the charged fermion before and after emitting the collinear 
photon, the function 0(zj) summarizes the phase space cuts. The integration over z% in Eq. (|5ip is 
constrained by the phase space cuts Q(zi) and cannot be performed analytically. Consequently the 
fermion mass singularities are not fully canceled in the combination of virtual and real corrections 
and thus become visible. If the photon is treated inclusively, the integration over Zi will not be 
constrained by any phase space cut, and thus can be performed analytically, leading to Eq. (|50p . 



B.2 Dipole subtraction 

In this approach [46,48,49] one constructs an auxiliary function which contains the same sin- 
gularities as the real bremsstrahlung integrand. Subtracting this auxiliary function from the 
bremsstrahlung integrand thus cancels all soft and collinear singularities and the difference can 
be integrated numerically, even in the singular region. In this numerical integration no regulators 
are needed for the soft and collinear singularities. The auxiliary function can then be integrated 
analytically (regulators required) and re-added to the original integral. Within the subtraction 
method there is no singular contribution involved in the numerical integration. Hence for compu- 
tations within this method, the statistical uncertainty is smaller than that of the slicing method, 
in which the singular contributions are present in the numerical integration. The auxiliary func- 
tion must possess the same asymptotic behavior as the original integrand in the soft and collinear 
limit, and has to be simple enough to be integrated over the singular regions analytically. In our 
case, soft and collinear singularities occur only in the final state. As the masses of the final state 
fermions can be consistently neglected, the expression of the auxiliary function is fairly simple 



|AW$4/ 7 )| 2 = - J2 Q* Qi 9ij\h, kj, k)\M born {^ iUj )\ 2 , (53) 



i,j=2 

#3 



where ki and k denote the respective momenta of final state fermions and photon, and the 
functions gfj b contain the soft and collinear singularities 



1 



2 



with the variables 



k%k k{kj /--v 

k^kj ~\~ k^k ~\~ kjk k^kj ~\~ kjk 

The mapping between the phase space of the radiative and non-radiative process, $4/ 7 and 
is defined as 

K = + W- -^-k] , k» = -i-fcj* (56) 
1 Uij 1 yij 

with all other momenta unchanged. The contribution of the auxiliary function should be com- 
puted analytically. After integrating over the photon momentum the result reads 

J ^4/ 7 |AW$4/ 7 )| 2 = -^Y,Qi Qi j d*Wi GtfiSi^Mbomi^ij)] 2 (57) 



with the function 



G!f(S ij )=L(S ij ,m*)-^ + l, (58) 



where L(Sij,mf) has been defined in Eq. 
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Figure 10: Results for the partial decay width of h — > WW* —> e~v e ^x + v^ in three different 
benchmark scenarios (the contribution due to the third generation fermion and sfermion loop 
corrections to the H°WW coupling is not included. TB denotes tan/3). The left column shows 
the corrected partial decay width, while the right column shows the relative corrections. 
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Figure 11: Contribution to the partial decay width of h° — > WW* — > e^v e ^ + v^ due to the 
correction to the H°WW coupling from the third generation fermions and sfermions in three 
different benchmark scenarios (TB denotes tan /3) . The left column shows the corrections to the 
partial decay width, while the right column shows their relative size. 
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Figure 12: Results for the partial decay width of h — > ZZ* — > e~e + fi + [i~ in three different 
benchmark scenarios (the contribution due to the third generation fermion and sfermion loop 
corrections to the H°ZZ coupling is not included. TB denotes tan /3). The left column shows the 
corrected partial decay width, while the right column shows the relative corrections. 



31 




Figure 13: Contribution to the partial decay width of h — > ZZ* — > e~e + due to the 
correction to the H®ZZ coupling from the third generation fermions and sfermions in three 
different benchmark scenarios (TB denotes tan/3). The left column shows the corrections to the 
partial decay width, while the right column shows their relative size. 
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Figure 14: Invariant mass distribution of in the decay h° — >■ e f e /^ + i/„ (no photon combina- 
tion, the contribution from the H°WW vertex correction is not included) in the m? 
with tan /3 = 30, M A o = 120 GeV (left) and tan /3 = 30, M A o = 400 GeV (right). 
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Figure 15: Relative correction to the invariant mass distribution of fi^u^ in the decay h — > 
e~v e n + v jjL (the contribution from the H°WW vertex correction is not included) in the m™ ax 
scenario, with tan/3 = 30, M A o = 120 GeV (left) and tan/3 = 30, M A o = 400 GeV (right), 
"com." and "no-com." indicate the results with and without photon combination, respectively. 
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Figure 16: Relative contribution to the invariant mass distribution of /x + z^ in the decay h c 



e v e /j, + i , f j i from the higher order final state radiation ("higher order FSR") and from the H°WW 

> imix scenario, with tan f3 



vertex correction ("3rd generation only") in the 
(left) and tan /3 = 30, M A o = 400 GeV (right). 
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Figure 17: Invariant mass distribution and relative correction to the invariant mass distribution 
of /U+f^ in the decay h° — > e~^ e /i + z/ At (no photon combination) in the no-mixing scenario (upper) 
and the small-a e fj scenario (lower), with tan/3 = 30, M A o = 400 GeV. 
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Figure 18: Invariant mass distribution of fi + fi in the decay h° — > e e + fi + fi (no photon combi 
nation, the contribution from the H°ZZ vertex correction is not included) in the 
with tan/3 = 30, M A o = 120 GeV (left) and tan/3 = 30, M A o = 400 GeV (right). 
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Figure 19: Relative correction to the invariant mass distribution of fi + fi in the decay h° — > 
e~e + n + [jL~ (the contribution from the H°ZZ vertex correction is not included) in the m™ ax 
scenario, with tan/3 = 30, M A o = 120 GeV (left) and tan/3 = 30, M A o = 400 GeV (right), 
"com." and "no-com." indicate the results with and without photon combination, respectively. 
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Figure 20: Relative contribution to the invariant mass distribution of /U + /x~ in the decay h — > 
e~e + [i + fi~ from the higher order final state radiation ("higher order FSR") and from the H°ZZ 
vertex correction ("3rd generation only") in the m™ ax scenario, with tan/3 = 30, M A o = 120 GeV 
(left) and tan /3 = 30, M A o = 400 GeV (right). 



-3l 

dAW 10 ' 1 S[%] 




Figure 21: Invariant mass distribution and relative correction to the invariant mass distribution of 
fi + n~ in the decay h° — > e~e + fi + n~ (no photon combination) in the no- mixing scenario (upper) 
and the small-a e ff scenario (lower), with tan/3 = 30, M A o = 400 GeV. 
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